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Abstract 

The paper develops the method for construction of families of par- 
ticular solutions to some classes of nonlinear Partial Differential Equa- 
tions (PDE). Method is based on the specific link between algebraic 
matrix equations and PDE. Admittable solutions involve arbitrary 
functions of either single or several variables. 

1 Introduction 

Analysis of nonlinear Partial Differential Equations (PDE) is severe problem 
in mathematical physics. Many different methods have been developed for 
analytical investigation of nonlinear PDE during last decades: Inverse Scat- 
tering Problem [1, 2, 3, 4, 5, 6, 7], Sato theory [8, 9, 10, 11], Hirota bilinear 
method [12, 13, 14, 15], Penleve method [16, 17, 18], <9-problem [19, 20, 21], 
with some generalizations [22, 23, 24, 25, 26], quasi-classical version of the 
(9-problem [27, 28, 29]. A wide class of PDE (so-called completely integrable 
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systems) has been studied better then others. Nevertheless, there are many 
methods which work in nonintegrable case as well: [12, 13, 14, 15, 16, 17, 18]. 

We represent the method for construction of the families of particular 
solutions to some classes of multidimensional nonlinear PDE. This method is 
based on general properties of linear algebraic matrix equations. Essentially 
we develop some ideas represented in the ref.[8] and recently in the ref. [11]. 

General algorithm is discussed in the Sec. 2. We consider systems, ad- 
mitting solutions depending on arbitrary functions of either single or several 
variables. Sec. 3 represents some examples among which are Kadomtsev- 
Petviashvili equation (KP) and Devi-Stewartson equation (DS). Conclusions 
are given in the Sec. 5. 



2 General results 

The algorithm represented in this section is based on the fundamental prop- 
erties of linear matrix algebraic equation, which is written in the following 
form: 



= (1) 

where \& = {tpij} is iV x N matrix, U and $ are N x M matrices. Let us 
recall these properties. 



1. If ^ is nondegenerate matrix, i.e. 

det ^ 0, (2) 

then equation (1) has unique solution which can be written in the 
following form: 

U = ^- l $. (3) 
Only nondegenerate matrices \I/ will be considered hereafter. 

2. (consequence of the previous property) If $ = and condition (2) is 
held, then the equation (1) has only the trivial solution 

U = 0. (4) 
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3. (superposition principle) Consider the set of K matrix equations with 
the same matrix ^: 



W i = $ il i = l,...,K. (5) 

Then for any set of scalars b k (k — 1, . . . , K), function U = J2k=i ^kUk 
is solution of the following matrix equation 

K 

W = Y,*k*k. (6) 

k=i 

4. (consequence of properties 2 and 3) If columns $j are linearly depen- 
dent, i.e there are scalars a k , k — 1, . . . , K, such that 

K 

J2 a k®k = 0, (7) 

fe=i 

then 

K 

Y^a k U k = 0. (8) 

k=l 

Note that analogous properties of linear integral equation have been used in 
the classical dressing method based on the <9-problem [19, 20, 21]. 

We will use two sets of variables, which will be introduced in the functions 
$ and \&: x — (x±, . . . , xq) , t = (ti, t 2 , . . . ), where Q = dim(x). The next 
statement follows from the above properties of the linear equations. 

If there is transformation T, which maps the nonhomogeneous equation 
(1) into the homogeneous 

W(U) = 0, (9) 

then 

U(U) = 0, (10) 

or 

$ = q?U ^ = W(U) =>• U(U) = (11) 
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If U depends on U and its derivatives, then equation (10) represents the 
matrix PDE for U. 

We will see that transformation T is not unique. One has the manifold 
of transformations T 

Tj G T : $ = W % = =>• 17,(17) = (12) 

which is uniquely defined by the equations introducing variables x and t in 
the matrices <3> and Namely, 

* Xk = *B k + QC k , k = l,...,Q (13) 

(B k and are N x N and M x N matrices respectively) and 

Q 

M^ = 0, M^ = 0, M2 = d ta +d a , J> fc = n, (14) 

fc=i 

where n is order of differential operator. We use Q-dimensional vector sub- 
script a = (ai, . . . ckq) and notation <9 a = n^Li^"/- Systems (13) and (14) 
should be compatible, which leads to relations among matrices Bj and Cj 
(see the next subsection). 

Each operator M™ defines transformation T a G T. In fact, let us apply 
M™ to both sides of the eq. (1) and use eqs. (14): 

= M™$ = (d ta V)U + Vd ta U + d a (WJ) = -(d a V)U + t>d ta U + (15) 
«i a Q / Q \ / Q \ / Q 

ii=0 jQ=0 \fe=l / \n=l / \m=l 

= 777 — ~T77 are binomial coefficients. Due to the eqs.(l) and (13) , one 
k\[n — k)\ 

has the following expression for d Xj ^>: 

d Xj * = *V e j, V e3 =B 3 + UC v e J = (0 1 __0, 1,(^^0). (16) 

j'-i Q-i 

Analogously we can write 

= VVp, (17) 
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where Vg can be found recursively. In fact, one has 

= n,,a, h , Jq) * = i; , , Q) ) = (is) 

(^^)%,..., /3j _i,..., /3q) + ^(^V^,...^.-!,...,^)) = 

i.e. 

V/3 = V^-V^,...^.-!,...,^) + d^V^,...,^-!,...,^), (19) 

^0 = ^(0,...,0) = T N,N- 

Q 

Owing to the eq.(17) we can replace factor ^IIn=i i R (^) with *&V(j 1>mmm j Q ), 

which results in 

= -(d a ^)U + ^d ta U+ (20) 

ai a Q / Q \ 

* e ■ • • e ( n c s ) %.-.*>> n ^ = 
= v + e ' (n c s) %--^) n ( 2i ) 

3i,—3Q \k=l / m=l 

where means the sum over all values ji, . . . , jg such that (X^m=i -7™) < 

3U—,3Q 

n The associated system of nonlinear PDE has the form (see (10)) 

^ + E [U C tj V (n,-,3 Q ) II ^~ jmU = ( 22 ) 

Ji,— 3Q \k=l J m=l 



2.1 Compatibility of the system (13) and (14) 

Studying the compatibility condition for the system (13) we concentrate on 
two particular cases 

1. compatibility condition produces the additional differential equation 
for the function U and Q — 1 differential equations for the function ip. 
In this case (f> (and U) involve N arbitrary scalar functions of single 
variable. 
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2. compatibility condition produces only the differential equation for the 
function U. In this case (f> (and U) involve N x M arbitrary scalar 
functions of Q variables. 

Intermediate case (i.e. p and U depend on arbitrary functions of P (1 < 
P < Q) variables) is also possible, but it is not regarded in this paper. 

2.1.1 Matrix U depends on N arbitrary scalar functions of single 
variable 

We start the analysis with splitting the eq. (13) into two equations by using 
the following structure of the matrices: 



B s = 








_ o N 


—M,M 


B 2j _ 


1 


Cq\ = 




3 = 


1,..., 


Q 



(23) 
(24) 



where <p is N x M, x is N x (N - M); bj, C 0j are M x M, B Xj and C Xj are 
M x (N — M); B 2j are (N — M) x (N — M) matrices for all j = 1, ... , Q; 
hereafter Oa,b with arbitrary A and B is Ax B matrix of zeros; Im is M x M 
identity matrix. Having matrices with given structure one can split the 
system (13) as follows: 

r V ■ «I>r 0; . (25) 

Xxj = xB 2j + ^Bij + Qdj (26) 
The eq.(25) with j = 1 defines $ in terms of ip (remember, that Coi = Im)' 

$ = ip Xl - ipb 1: (27) 

Substituting this expression into the eq.(25) with j > 1 one gets the overde- 
termined system of equations for ip : 

<p Xj = VxiCoj + ip{bj - hCoj), j > 1. (28) 

Compatibility of this system with different j leads to the following relations 
among matrices: 

[C ok ,C Oj ] = 0, (29) 
[bj, Cok] — [bk, Coj] + CokbiCoj — CojbiCot = 0, 
[bj - hCoj, b k - hCok] = 0. 
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From another point of view, eqs. (1), (13) and (25) lead to the additional 
differential equation for U: 

E M [bi, bj) + UCoibj + U Xi C 0j + BiUCoj + UCiUC 0j = (30) 
U C jbi + U Xj C 0i + BjUCoi + U CjUC 0i , 

where Em is matrix,composed of first M column of the iV-dimensional iden- 
tity matrix I N . So, eq. (25) describes one more transformation T com P e T. 

The system (26) is overdetermined system defining x m terms of (p. Its 
compatibility condition leads to the following relations 

A 3 x + A 2 <p xixi + A x ip xi + A ip = 0, (31) 
A 3 = [B 2i ,B 2j }=0 
A 2 = CoiCij — CojCu = 0, 
A n = A n (Coi, Coj, Cu, Cij, Bj) = 0, n — 0,1. 

We don't represent expressions for A and because they are too compli- 
cated. But in the particular case, when bj = 0m,m and C\j = 0m,n-m the 
above system has simpler form 

[B 2i ,B 2j ] = (32) 
CoiCij — CojCu = 0, 
CuB 2 j + CoiBij — C\jB 2 j — CojBu = 
BuB 2 j — B\jB 2i = 0. 



To satisfy both systems (13) and (14) ip should be written in the form: 



<P 



c(ki) exp 



Q 



i=i 



dki, 



(33) 



where c(fci) (N x M matrix function of argument) and kj satisfy the disper- 
sion relation associated with eq.(28), u a satisfy the dispersion relation for 
eq.(14). One can see that the above expression for <p> (and U due to the eq. 
(3)) involves TV arbitrary functions of single variable. Nonlinear equations 
generated by the operators M™ have the form (22) 
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2.1.2 Matrix U depending on N x M arbitrary scalar functions of 
Q variables 

In this section we consider matrix equation (1) in the form 

$ = q?U, $ = [ $1 • • • $ P ], U=[ XJx ■■■ Up], or (34) 

$ fc = W fc , k = l,...,P, (35) 

where <£>j and Uj are N x M, is N x N matrices. Let matrix Ck be of the 
form: 



®M(k-l),N 

Ck 

®N-Mk,N 



(36) 



where Ck is M x N matrix. So that eq.(13) can be written as follows: 

* Xfc = *B k + $ k C k , k=l,...,Q, (37) 

variables t n are introduced by the same eq.(14). Assume the following struc- 
ture of matrices: 



V=[<p\x], C k =[l M \C lk ], B k 



b k 


B\k 


®M,N-M 


B 2 k 



(38) 



where b k , are M x M matrices, C\ k and B\ k are M x (N — M) matrices, B 2k 
are (N — M) x (N — M) matrices for all k — 1, . . . , Q. Now we can split the 
eq.(37) into the following pair of equations 



<p Xk = (pb k + $ fe , 
Xx k = xB 2k + <pB lk + $ fc Cifc. 



(39) 
(40) 



Equation (39) defines function in terms of tp. At the same time, <& k can 
be eliminated from the eq.(39) due to the eq.(35), so that eq.(39) represents 
overdetermined system for p with compatibility condition (compare with 
eq.(30)) 



E M [bi, bj] + U J Xi + U% + B0 + \PC0 = 

ui. + u'bt + Bjir + wCjir, 



(41) 
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So, additional transformation X com P G T is based on the system (39). 

Compatibility of the system (40) gives relations among matrices Bj and 
Cj : for all j and k 



Cij — Cik = 0, 
[B 2 j,B 2 k] = 0, 
Bij — bjC\j + C\jB2j = 0. 

Expression for (p satisfying both systems (13) and (14) is following: 



(42) 



c(ki,. . . ,k Q )exp 



i=i 



dki . . . dkg, 



(43) 



where c is N x M arbitrary matrix function of arguments; u a satisfy the 
dispersion relations for (14). Due to the formula (3) and (43), U depends on 
N x M arbitrary functions of Q variables. Nonlinear equations generated by 
the operators M™ keep the form (22) 



3 Examples 

We consider two hierarchies of PDE, admitting solutions depending on iV 
functions of single variables (see subsection 2.1.1). Inside of them are classical 
KP and DS hierarchies. It is important to note that iV can be arbitrary 
integer for later, which is not true in general. 

3.1 KP hierarchy 

Let M = 1, Q = 1, M n = d tn + d£, n = 2, 3, . . . . Operators M n generate the 
following hierarchy (see eq.(22)): 

n 

U tn + J2c:V n ^dlU = 0, (44) 
V = I N ,;V 1 = B + UC, V n = KK-i + c^K-i, 

where B = B\ and C = C\. In this case one has only one eq.(13) and no 
compatibility condition for it, i.e. matrices B and C are arbitrary matrices 
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having structure (23) and (24). If 








1 





... 










1 


... 


B = 











1 ... 



c= [ 1 1 ... ] , u 



U N 



, (45) 



iV is arbitrary, then system (44) represent KP hierarchy [8 
written for the function u = U\ x using eq.(44) with n = 2, 3: 



KP can be 



1 3 

U t + ^UXXX + ^d~ l Uyy + 3UU X = 



(46) 



3.2 DS hierarchy 

Let M = 2, Q = 2, x\ = x, x<i = y. Consider operators M n = d tn + d™. This 
operators generate the same hierarchy (44) with appropriate dimensions for 
matrix U, 

V 1 = B 1 + UC 1 , (47) 

and additional nonlinear equation (30): 

U X C 02 + B 1 UC 02 + UC\UC 02 = U y C m + B 2 UC 01 + UC 2 UC 01 (48) 

with relations among matrices Bj and Cj given by (29) and (31) (or (32)). 
DS hierarchy corresponds to the following choice of matrices Bj and Cf 



B 1 









1 





... 











1 


... 














1 ... 



Ci 














1 





... 










1 








... 


B 2 = 

















1 ... 
















1 


... 



" 1 





... 


, c 2 = 


' 


1 


... 





1 


... 


1 





... 



u 



Ui Vi 



UN V r , 



(49) 
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N is arbitrary. It is simply to check that conditions (29,32) are satisfied. 
Let us consider the eq.(44) with n = 2 together with eq.(48). The following 
change of variables 

ui = |(r - q + 2w x ), u 2 = ^ (r + q + 2w y ) , (50) 

1. 1. 
«i = -^{-r - q + 2wy), v 2 = -(-r + g + 2w x ). 

transforms them into 

r t 2 - r xy - 2rw xy = 0, q t2 + q xy + 2qw xy = 0, w xx - w yy = qr, (51) 

which results in DS after reduction r — ip, q — ip, t 2 — it, where i 2 = — 1, 
bar means complex conjugated value 

4 Conclusions 

The suggested version of the dressing method is one more form of represen- 
tation of the systems of nonlinear PDE, which admit an infinite number of 
commuting flows, generated by the operators M™. Among them are classical 
completely integrable systems of equations. The feature of these system is 
that their solutions admit any number of arbitrary scalar functions of single 
variable (N is arbitrary for them). One can see that there is another type 
of systems, which admit an infinite number of commuting flows, while iV is 
given. We don't know interesting examples for this case. We also have not 
found the systems, whose solutions would admit any number of arbitrary 
scalar functions of several variables: for all nontrivial examples this number 
has to be fixed. These aspects will be studied later. 

The work is supported by RFBR grants 01-01-00929 and 00-15-96007. 
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